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Abstract 

We present type IIB supergravity solutions wliicli are expected to be dual to certain 
Lifshitz-like fixed points with anisotropic scale invariance. They are expected to describe 
a class of D3-D7 systems and their finite temperature generalizations are straightforward. 
We show that there exist solutions that interpolate between these anisotropic solutions in 
the IR and the standard AdS^ solutions in the UV. This predicts anisotropic RG flows 
from familiar isotropic fixed points to anisotropic ones. In our case, these RG flows are 
triggered by a non-zero theta-angle in Yang-Mills theories that linearly depends on one of the 
spatial coordinates. We study the perturbations around these backgrounds and discuss the 
possibility of instability. We also holographically compute their thermal entropies, viscosities, 
and entanglement entropies. 
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1 Introduction 



The AdS/CFT correspondence has provided a very powerful and successful paradigm to analyze 
relativistic and isotropic fixed points in various quantum field theories [H [21 El S] . On the field 
theory side, they are described by {d+ l)-dimensional conformal field theories and are invariant 
under the homogeneous scaling transformation {t, xi, X2, ■ ■ ■ , x^) — > (At, Xxi, Ax2, • • • , Ax^). On 
the gravity side, they are equivalently described by gravity on a (d + 2)-dimensional AdS space 

ds^ = [-dt" + Ildx^^ + TT- (1-1) 

It is natural to try to extend the AdS/CFT correspondence to a holography for the following 
anisotropic spacetime 

ds^ = Ldt' + ^2 d^'] + E + (1.2) 

V i=l J j=P+l 

where < p < d — 1, and the parameter 1) measures the degree of Lorentz symmetry 
violation and anisotropy. Since the metric (jl.2p is invariant under the scaling (t,Xi,yj,r) — > 
{X^t, X^Xi, Xyj, j), we expect that on the field theory side it is dual to a fixed point which is 
invariant under the scaling transformation 

{t,Xi,yj) {XH,X'xi,Xyj). (1.3) 

Notice also that by a coordinate redefinition = p, we can rewrite the metric (|1.2p into another 
illuminating form (after rescaling {t,Xi,yi) accordingly) 

ds^ = (-dt^ + dxi\ + pi E + (1-4) 

V i=l J i=P+l ^ 

Thus we can equally argue that the dual background is invariant under an anisotropic scaling 
transformation {t,Xi,yj, p) {Xt, Xxi, X^yj, j), where the yj directions are responsible for the 
Lorentz symmetry violation and anisotropy. 

In general, fixed points with the anisotropic scaling property (|1.3p are called Lifshitz(-like) 
fixed point^ (see e.g. the textbook [7] for a brief review). This generalization of AdS/CFT 
correspondence to Lifshitz-like fixed points ()1.2p was first proposed and analyzed by Kachru, 
Liu and Mulligan [8j in the particular case of p = 0. The simplest case with p = represents 
non-relativistic fixed points with dynamical critical exponent z, which appear in many examples 



* The most standard example is the free scalar field theory with z = 2, known as the Lifshitz model. Anisotropic 
fixed points in interacting field theories in general can have z ^ 2 [5||6]. Even though the original Lifshitz fixed 
points were found in anisotropic magnets where three critical lines meet, in this paper we simply define Lifshitz-like 
fixed points as any fixed points which have anisotropic scale invariance. 
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of quantum criticality in condensed matter physics (see the references in [8j ) . See [21 [THl IH [121 
[T3l [H] for further progress on holographic aspects of this topicj^ 

The other cases where 1 < p < d — 1 are not only generalizations of p = case but can 
also be interpreted as space-like anisotropic fixed points (see also [13j) as is clear from the 
expression (|1.4|) . Lifshitz fixed points with space-like anisotropic scale invariance appear in 
realistic magnets such as MnP and the axial next-nearest-neighbor Ising model [5]. They are 
also realized in models of directed percolation [7]. 

To understand holographic duals of such gravity backgrounds, it is the best to embed them 
into string theory, where microscopic interpretations are often possible by using D-branes. How- 
ever, so far there has been no known embedding of (jl.2p in string theory. Motivated by this 
circumstance, in this paper we will construct such anisotropically scaling solutions in type IIB su- 
pergravity. We mainly focus on the backgrounds generated by intersections of D3 and D7 branes. 
They correspond to the choice p = 2 and d = 3 and are expected to be non-super symmetric. 
This restriction is imposed not only for the tractability of the supergravity analysis, but is also 
due to another motivation, namely to construct back-reacted D3-D7 solutions that are dual to 
the pure Chern-Simons gauge theory in the second setup of |22jJ. In the end, we find a class 
of solutions with the exponent z = 3/2. We also extend them to black brane solutions dual to 
finite temperature theories. 

Furthermore, we show that there exist solutions which interpolate between our anisotropic 
solutions and the familiar AdS^ x solutions. We also construct their numerical solutions. 
The holography suggests that our Lifshitz-like fixed points can be obtained from various four- 
dimensional CFTs including = 4 super Yang-Mills via RG flowjfl. These flows are triggered by 
the relevant and anisotropic perturbation which gives a non-zero 9 parameter (i.e. the coefficient 
in front of the topological Yang-Mills coupling F A F) that depends linearly on one of the three 
spatial coordinates i.e. 6 oc x^. Notice that when X3 is compactified, the perturbation induces 
the Chern-Simons coupling f A A F + ^A^ as in [22], which becomes relevant in the IR. 

This paper is organized as follows. In section [2l we present solutions dual to a class of 
Lifshitz-like fixed points based on D3-D7 systems with their black brane generalizations. In 
section [31 we show there exist interpolating solutions which approach the Lifshitz-like scaling 
solutions in the IR and the standard AdS^ solutions in the UV. In section [H we holographically 
calculate the shear and bulk viscosity. In section \5\ we compute their holographic entanglement 
entropies and discuss how the scaling behaviors of the entanglement entropies depend on the 
direction along which the sub-systems are delineated. In section [HI we study the perturbations 
around these backgrounds and discuss the instabilities. In section [71 we present anisotropic 
solutions based on D4-D6 systems. In section [HI we summarize our conclusions. 

^ Gravity duals of another types of fixed points with non-relativistic scaling symmetry [15] have also been 
studied especially for systems with non-relativistic conformal invariance [161 1171 1181 1191 1201 121] . 

^Here the Lifshitz-like fixed points are realized in the IR limit. It is also intriguing to consider opposite RG 
fiows, where IR fixed points become relativistic and isotropic z = 1 as in [51 110). 
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2 Holographic Duals of Lifshitz-like Fixed Points in Type IIB 
String 

In this section we will present the main result of this paper. We will construct new solutions 
in type IIB supergravity with RR 5-form and 1-form fluxes whose Einstein metrics enjoy a 
nice scaling property. Since their scaling is anisotropic as opposed to the well-known AdS^ 
background, we argue that they are dual to Lifshitz-like fixed points described by certain D3-D7 
systems. 

2.1 Type IIB Supergravity 

The IIB supergravity action Sjjb = f £ in the string frame is defined by the Lagrangian 
(we follow the convention in j23j ) 

g-20 1 1 

£ = V^e-^'^{R + 4dM(l)d^(l>) - -^H3 A *Hs - -F, A *Fi - -F3 A *F3 

-^Fs A*F5-^C4A/73AF3, (2.1) 

where Fi = dx, F3 = F3 - XH3, and F5 = F5 - A H3 + \B2 A F3. We set a' = I therefore 
2Kfo = (27r)7. 

The fluxes obey the equations of motion: 

d*Fi = *F3AH3, d* F3 = -i?3 A F5, d*F5=H3AF3, 
d(e"2'^*i?3) = Fi A*F3-hF3 AF5. (2.2) 

plus the Bianchi identities: 

dH3 = 0, dFi = 0, dF3 = H3AFi, dF^ = H3AF3, (2.3) 

and the self-dual constraint for F^: 

*h = h- (2.4) 

The dilaton equation of motion is 

R + 4VmV^V - ^HmnpH^'^'' - 4Vm<^V^V = 0. (2.5) 
And the Einstein equation becomes 

Rmn + '2'^m'^n4> + ^gMNA 

= ^^HmabHj,^^ + ^e^'^FuF^ + -^c^'^FmabF^^^ + ^e't'FMABCDFN^'''''' , 

(2.6) 

where 

A ^ e'^duxd^'x + ^e^^FABcF^^^ + ^e^'^F^^ci.EF^^^^^. (2.7) 
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2.2 D3-D7 Ansatz 



We start with the (string frame) metric ansatz that preserves the three-dimensional Lorentz 
symmetry S0{2, 1): 

ds^ = e^'('-\-dt' + dx^ + dy^) + e2'^M+2'^(^)d«;2 + e^^^^-^^^^Ur^ + e^^^'^V^d^^. (2.8) 

We require the five-dimensional compact manifold X5 to be a unit-radius Einstein manifold with 
the same Ricci curvature as the unit-radius S^, i.e. it satisfies 

Rap = ^9a(3- (2.9) 

The simplest example of X5 is obviously the unit radius sphere 5^. The self-dual 5-form and 
1-from fluxes are given in terms of constants a and (3 by 

F5 = a(17x5 + *^^X5), (2.10) 
Fi = dx = (3dw, (2.11) 

where x is the axion field (i.e. the RR 0-form potential) and f^Xs is the volume form of X5. The 
fluxes (j2.10p and (j2.1ip satisfy the equations of motion (|2.2p . We also assume that the dilaton (p 
only depends on r and both 3-form fluxes (i/3 and F3) vanish. Our ansatz, which looks rather 
different from [8j , is motivated in part by an attempt to construct back-reacted solutions of the 
D3-D7 intersecting systems introduced in [22j, as will be explained in detail later. 

Under this ansatz, the equations of motion for the metric and the dilaton ( (j2.6p and ()2.5p ) 
are summarized as follows: 

02 2 

[(a + hye^']' = _^e-2a-h+3fe+6c^5 ^ ^^-4c+3b+h^-5^ 

[(c + logr)'e2^]' = 4e''"'" + ^g-^'^-'^+^^+^'^rS - ^g-^^+^^+V^^^ 



f2 4 

/ 

"4 " '4 



2z" + c" - a" + 2{z'f + ^{h'f + a'h' + 2{c'f + (- + a')c' + ^{b'f - + ttI = «• 

2 r 2 2r^ 

Here we have defined 

^ = lb + l^°Sr + a + 2c+^h-(j). (2.12) 

The derivative of a function / with respect to r is denoted by f'{r). An observation, which will 
be useful in the next section, is that a linear combination of the first four equations gives 

[(26 - 2a - - 2/i)'e2^]' = 0. (2.13) 
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2.3 D3-D7 Scaling Solutions: Holographic Duals of Lifshitz-like Fixed Points 

Since we are looking for scaling solutions (namely solutions invariant under scale transforma- 
tions), we require all metric components in ()2.8p to be power functions of r. In other words, 
the functions a, b, c, z and 4> are all logarithmic functions of r. For such a scaling ansatz, the 
equations of motion (j2.5p and (|2.6p reduce to algebraic equations and the solution is easily found 
to be: 



a(r) 
h{r) 



1 12 7^ 

-log— -log ^s, b{r) = log r + bo, c(r) 

2 11 



-1 + ^ ) logr + Co, 



56 



log r + log 6 + ho, 



(r) = ^ log r + 



ho-CQ- 



(2.14) 



where 6o)Co,/io,0o and ^ are arbitrary constants. ^ corresponds to the degrees of freedom 
of the reparameterization of r, while 6o and ho correspond to the rescaling of the {t,x,y,w) 
directions. 

Without loss of generality, we choose 



6 = 1, ^0 
and the solution in the string frame reads: 



1 , 11 

CO + 2 log 



hn 



, 11 

Co + log— , 



dr 



rs{-dt + dx + dy ) + r^dw H 



r 3 



+ Rlrzds\^, 



where i?^ = jji?^ = 6^"^°. And the dilaton scales with r as 



r 3 e 



4>o 



(2.15) 



(2.16) 



(2.17) 



where e"^o = 

Since the dilaton depends on r non-trivially, it is helpful to discuss the metric in the Einstein 
frame. Indeed, later we will see explicitly that a large class of scalar fluctuations around this 
solution can be described by Klein-Gordon equations on curved spacetimes based on the Einstein 
frame metric instead of on the string frame metric. The above solution in the Einstein frame is 



ds\ = 



r'^{-df + dx"^ + dy"^) + r^^dw^ + 



dr^ 



+ R^dsV 



where the radii 



i?2 = 

11 



■^0 

e 2 



-2co 



The metric ()2.18p is invariant under the scaling 



2 T 

{t, X, y, w,r) ^ [ At, Ax, Xy, X~iw, - 



(2.18) 
(2.19) 

(2.20) 
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and therefore is expected to be holographically dual to Lifshitz-like fixed points with space-like 
anisotropic scale invariance. Note that the metric ()2.18p is equivalent to ()1.4p with z = 3/2, 
p = 2 and d = 2>. 

2 

By redefining the radius coordinate p = rs and rescaling {t, x, y, w) accordingly, we can 
rewrite the metric (|2.18p into another illuminating form 



p^i-dt"^ + dx"^ + dy"^) + p'^dw'^ + 



dp2 



(2.21) 



This can be regarded as gravity duals of Lifshitz-like fixed points with z = 3/2. It coincides 
with the metric (|1.2p with p = 2 and d = 3. 



2.4 Holographic Interpretation in terms of D3-D7 System 

Since our solution (j2.18p is sourced by the RR 5-form (j2.10p and 1-form flux (j2.1ip . we expect 
it to be interpreted as a D3-D7 system in string theory. When we compactify the w direction 
such that w ^ w + L and place D3-brane along the (t, x, y, w) directions and k D7-branes 
along the {t,x,y, X^) directions: 



M4xS^ X X5 


t 


Xl 


y r 


w 


Si S2 S3 S4 S5 


N D3 


X 


X 


X 


X 




k D7 


X 


X 


X 




X X X X X 



these N D3 and k D7 branes can source the desired RR 5-form and 1-form fluxes with 

{27r)^N 



a 



Vol,A,)- P-^^' 

This brane configuration is the same as the one constructed to model the fractional quantum 
Hah effect in [22]. 

The number of the D3-branes determines the radii R and R in the scaling solution (j2.18p : 



Ri' 



12 ~. 



11 



-R' 



-N. 



(2.23) 



V Vol(X5)" 

For X5 = (whose volume is vr^), R^ = jf-R^ = 2^-kN . The number of the D7-branes gives 
the string coupling at r = 1: 

e^o = (2.24) 

Now we say a few words about the field theory living on this D3-D7 system. We take 
X5 = to simplify the arguments. If we start with N D3-branes, whose low energy theory is 
the four-dimensional = 4, SU{N) super Yang-Mills theory, then the additional k D7-branes 
will source a non-trivial axion field x = x^' '^^^ich in turn induces a tn-dependent 9 term (i.e. 
the topological term) of the Yang-Mills theory 

^ j x(u^)TrF A F. (2.25) 



7 



For finite [3 and w-direction is compactified. After integrating over w, tlie 4D topological 
term (|2.25p becomes a 3D Chern-Simons term at level k: 



k_ 
in 



Tr 



2 o 

aaf + -a^ 

3 



(2.26) 



Now we have two choices of the boundary condition for the u;-circle: periodic or anti-periodic. 
If we impose the anti-periodic one, all fermions will become massive. This breaks all supersym- 
metries and gives masses to scalar fields through quantum corrections. In the IR limit, only 
a pure Yang-Mills term is left of the original 4D = 4 super Yang-Mills part of the action. 
Since in the IR limit, the Chern-Simons term dominates this Yang-Mills term, the final three- 
dimensional theory is a pure Chern-Simons theory. In [22j, this D3-D7 system was constructed 
to holographically model the FQHE precisely because it flows to the pure Chern-Simons gauge 
theory in the IR. In this model, the AdS/CFT correspondence in the IR limit manifests itself 
as the level-rank duality of the pure Chern-Simons gauge theory. 

On the other hand, if we take /c — > oo (and simultaneously L — > oo) while keeping /? finite, the 
tD-direction is non-compact and the field theory is four-dimensional. Even though the interaction 
(j2.25p looks non-local at the first sight, its contribution to the equations of motion is actually 
local. This remarkable property occurs only when x(^) is a linear function of w (as is the case 
here). 

One might still doubt any relations of our new background ()2.18p to the = 4 super Yang- 
Mills theory as it is not asymptotically AdS^. One might also worry that the dilaton (j2.17p blows 
up near the boundary r oo. However, as we will show in the next section, we can in fact 
construct solutions which interpolate between the AdS^ and our scaling solution (j2.18p . This 
interpolating solution can be considered as the dual of the RG flow between the two systems. 
Notice that this caps off the strongly coupled region of the scaling solution. We will also present 
anisotropic solutions for analogous D4-D6 systems in section [71 

2.5 Black Brane Solutions and Entropy 

One more interesting fact about our scaling solutions is that we can straightforwardly generalize 
them to black brane solutions which have regular event horizons. The metric in the Einstein 
frame is 



dsl = 



dv 

r\-F{r)dt^ + dx^ + dy^) + r^dw^ + 



+ R^ds\^, (2.27) 



where 



F{r) = 1-At. (2.28) 

T 3 

The constant represents the mass parameter of the black brane. The dilaton and RR fields 
remain the same. 
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Requiring the smoothness of the Euchdean geometry of (|2.27|) gives the Hawking temperature 

Th = ^,^^. (2.29) 

The Bekenstein-Hawing entropy is then 
where 7 is a numerical factor 

7 = 2! • 3^ • • vrt ~ 3.729 (2.31) 

and V2 represents the area in the (x, y) direction. The entropy (j2.30p is proportional to N'^ and 
thus is consistent with the planar limit of a certain gauge theory. 

Notice that the power 8/3 of temperature in (|2.30p can also be obtained from a simple 
dimension counting. From the metric (j2.27p . the coordinate w has the fractional dimension 2/3, 
while each of (t, x, y) carries the unit dimension. 



3 RG Flow in AdS5/CFT4 and Scaling Solution 

In the previous section, we find a new scaling solution of the D3-D7 system in type IIB su- 
pergravity. To clarify its physical interpretation, we will show below that we can construct 
interpolating solutions that approach the AdS^ x solutions in the r — > 00 limit (i.e. UV 
limit of the holographic duals) and the scaling solutions in the opposite limit r ^ 0. Then 
via the AdS/CFT correspondence, we can argue that the system dual to our scaling solution is 
connected to the one dual to the ^4^55 through the RG flow. 

3.1 Further Reduction of Equations of Motion 

To find the interpolating solution, we start with the general form (|2.8p . To simplify the problem 
we impose some extra constraints which are consistent with both the AdS^ and the scaling 
solutions. 

First, we can make the function a{r) vanish by a reparametrization of r: 

a{r) = 0. (3.1) 

Secondly, recall that we showed [(26 — 2a — — 2/i)'e^^]' = for generic solutions. In fact, both 
the AdS^ and the scaling solution satisfy a much stronger condition 

(26 - 2a - - 2h)' = 0. (3.2) 
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Since we are looking for a solution that interpolates between the AdS^ and the scaling solution, 
it is reasonable to impose p.2p as a simplifying ansatz, namely 

h{r)=b{r)-a{r)-^cPir) + ho, (3.3) 

where ho is a constant. Similarly, since both the AdS^ and the scaling solution have the nice 
property that in the Einstein frame the radius of is a constant, we will also impose this 
condition on our interpolating solution, namely we require 

(j){r) = 4c{r) + 4logr + 4>o, (3.4) 

where (po is a constant. These relations p.lh . (13. Sp and (|3.4p are the constraints mentioned and 
are assumed throughout this section. 

Under this ansatz, the Einstein frame metric becomes 

^2b{r)-2c{r)-fy _ ^2b(r)-6c(r) t t 

dsl = ' J \ -dt'+dx'+dy')+e-lf">^'^o ' dw'+e-'^^+e-'-^dsj,^, (3.5) 

and the equations of motion are greatly simplified: 

a = Ae-^\ (3.6) 
b" = ^ + -b' - 106'2 - 1^ + 246'c' - 8c'2, (3.7) 

c" = 4 + —5' - 66'2 _ ^ + l46'c' - 2c'2, (3.8) 

We can confirm that the derivative of the r.h.s of (j3.9p is vanishing if (j3.7p . (j3.8p and (j3.9p are 
satisfied. This means that the constraint (13.91) is consistent with (13.71) and (l3i 



3.2 Interpolating Solution between AdS^ and D3-D7 Scaling Solution 

Now the problem amounts to solving the system of two coupled first-ordered nonlinear ODEs 
()3.7p and (j3.8p under the constraint (|3.9p . First, notice that (|3.7p and (|3.8p involve only the 
derivatives of b and c, thus once we find a solution to them, we can simply choose the integration 
constants of b{r) and c(r) such that they satisfy the constraint ()3.9p — as long as it allows the 
r.h.s of (j3.9p to be positivejll Therefore essentially we only need to solve p.7p and (j3.8p . 

Next we redefine the radial coordinate r and the derivatives of the functions 6(r) and c(r) 
as follows 

,.logr. B(s,^§^, CW.|*i. (3.10) 

o log r a log r 



^However as we will show later this requirement is actually automatically satisfied by the interpolating solution 
we are looking for; thus it does not impose any additional constraint. 
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Then the equations (|3.7|) and (|3.8|) are simply a pair of first-ordered non-hnear ODEs: 

B = 2 + 2AB-l6C- lOB^ + 24.BC - SC^, (3.11) 

C = 4 + UB-AC-6B'^ + UBC-2C'^. (3.12) 

where B = A physical solution also needs to satisfy 

(95 - 8C - 8)2 - (33^2 + 48) < 0, (3.13) 



due to the constraint (13. 9p . 

The dynamical system (13.1ip and (j3.12p has four fixed points {B,C)* which can be classified 
into two pairs: 

(S,C)* = (±1,-1) and {B,Cf = {±^,±-^-l\ . (3.14) 



33' \/33 

Inside each pair, the two fixed points are related by a coordinate redefinition r — > ^ thus are 
equivalent. The fixed point 

(^,C')Ad55 = (l'-l) (3-15) 

corresponds the standard AdS^^ x solution^] 

Our scaling D3-D7 solution (at zero temperature) corresponds to the fixed point 

(S, C):,,ung = (4|, - l) (1.2185, -0.8259). (3.16) 



One can easily see that the metric ()3.5p with ()3.16p is equivalent to (j2.18p via the redefinition 
of radial coordinate r — > 

^^33/6^ Since the two fixed points with "— " sign are equivalent to the 
two with "+" sign and are disconnected from them, we will not consider those any further. 

Now let's study the behavior of this dynamical system ()3.1ip and ()3.12p . Near the AdS^ 
fixed point (j3.15p . the eigensystem of the linear perturbations (defined by B{s) = B*^^g^ + eb{s) 
and C{s) = C\^g^ + ec{s) for e^, Ec <C 1) is 

€6 = -20eb + 24ee, = -12e6 + 14e„ (3.17) 

and both eigenvalues are negative: Ai = —4 and A2 = —2; therefore the AdS^ solution is a stable 
fixed point as the system flows to the UV (i.e. r 00). 

On the other hand, the eigensystem of the linear perturbation near the D3-D7 scaling fixed 
point (defined by B{s) = B*^^^^^ + 7jb{s) and C(s) = C*^aiing + Vc{s) for T]b,Vc <C 1) is 

116 152 70 94 , , 



*The fixed point {B,C)* — (—1,-1) can be considered as its conjugate since they are connected under a 
coordinate redefinition r -. 
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In contrast to the stable AdS^ fixed point, this scaling solution fixed point has one negative 
eigenvalue (Ai = - v^+v^ ) aj^^^ Qj^g positive one (A2 = ^^-'•^ ). Therefore the fixed point 
corresponding to the scaling solution is unstable. Near the neighborhood of this scaling fixed 
point, there exist two special trajectories: one corresponding to the negative eigenvalue Ai and 
one to the positive A2. The fixed point behaves like a UV (resp. IR) fixed point when approached 
along the trajectory corresponding to the negative (resp. positive) eigenvalue. When the fixed 
point is approached along a generic direction, the trajectory only passes near it and then turns 
to flow to infinity — only one fine-tuned trajectory can reach the fixed point. 



Figure [T] shows the global behavior of this dynamical system ()3.1ip and (I3.12|) . Figure 1(a) 



is generated by numerical computations and the salient features are schematically highlighted in 
the hand-rendered Figure 1(b) The arrows point in the direction from the UV (s = 00) to the 
IR (s = —00). It is easy to see from the direction fields that the AdS^ fixed point (the green dot 
at (1, —1)) is a stable UV fixed point while the scaling fixed point (the blue dot at (—£5, ^5 — 1)) 

V 00 V 00 

is unstable. Recall that a physical solution needs to satisfy (|3.13p . The allowed region in the 
fiow diagram are between two hyperbolic lines given by 95 - 8C - 8 = ±^33^2 + 48. The black 
curve in Figure 1(a) is the one with the sign; and the other one with the "— " sign is its 
mirror image in the upper- left corner but is out of the range of Figure 1(a) It is clear that both 
the ^^55 fixed point and the scaling solution fixed point are in the allowed region. And since 
the AdS^ is a stable UV fixed point and there is no critical surface separating it from the scaling 
solution fixed point, there exists a trajectory emanating from the AdS^ fixed point and flowing 
to the scaling solution fixed point. Namely there exists a solution that interpolates between the 
AdS^ X X5 solution in the UV (r — > 00) and the D3-D7 scaling solution in the IR (r = 0). 

Now to solve the interpolating solution, we first choose the integration constant. First, (j)o 
is determined only by the 5-form flux: 0o = ~ log f = — log ( y^i^ -^^ ^ ^) ■ Then without loss of 

generality, we set Jiq = ^ and choose the boundary condition for (6, c) to be 

6(r)— >logr, c(r)—>— log r at r ^ +00, (3.19) 

b{r)^ -J= log r + bo, c(r)^ T-^-lj log r + &o at r ^ 0, (3.20) 



where 60 = | log ( Vol(x ) ~ H 1 ~ I Z^- ^hen the dynamical system (f3lT]l and (f3T^ 
plus the constraint 

- [{9B - 8C - 8)2 - (33^2 + 48)] = — ^ g-^M'^^+^^^^'^Vi^ < 0, (3.21) 

«(to&) 

and the boundary conditions (j3.19p determine an interpolating solution that approaches the 
^4^5*5 X solution 



dsl = 



r'^{-dt^ + dx^ + dy^ + dw^) H ^ 



+ R'ds\^, 



(3.22) 
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^^^^^^^^^^^^^ 
y- y- y- y ^ y^ 
^ ^ ^ ^ ^ y^^'''^/^/^ y 
^ ^ ^ ^ ^ ^y^^^^ ^ ^ 

^ ^ y ^ ^ ^ y 
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y-yyy'y'yyyyy'y' 
yyyyyyyyy- 
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(a) 



(b) 



Figure 1: The {B,C) flow diagram. The horizontal and vertical axes are B = b and C = c, 
respectively. The arrows point in the direction from the UV (r = oo) to the IR (r = 0). In the 
left figure, the blue dot at (—£5, ~ (1-22, —0.83) is the unstable fixed point corresponding 

to the scaling solution; the green dot at (1, —1) is the stable UV fixed point corresponding to 
the AdS^ solution. The green line running through the scaling fixed point corresponds to the 
negative eigenvalue Ai while the red line corresponds to the positive one A2. The black curve is 
given by 9B — 8C — 8 = \/33-B^ + 48 and corresponds to a pure D3 solution; the allowed D3-D7 
solutions are above this curve. 
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in the UV and the D3-D7 scaling solution 



dSE 



12 (]jT^ ^ 



(3.23) 
(3.24) 



in the IR. Here = e-<^o/2 ^ 2j^j^^^N and = e'^^^'R^. And the fluxes are 



throughout the system. 

We can easily solve (6, c) numerically for arbitrary fluxes. Systems with difl^erent flux numbers 
difi'er only in their speeds in approaching the fixed points. Figure [2] shows the system with 



/32 



64{ 



Vol{X5) 



)2Ar2 



1 as an example. Figure 2(a) graphs the behavior of (6, c), and more directly, 



Figure 2(b) presents the scalings of the (t, x, y)-directions, w-direction, and e'^ in the Einstein 



frame as the interpolating solution flows from the D3-D7 scaling solution in the IR (r = 0) to 
the AdS^ in the UV (r ^ oo). 



I (t,x,y) w e'^{phi}] 



-5 -4 -3 -2 



2 3 4 5 

logW 




(a) (B, C) of the interpolating solution. 
Figure 2: An interpolating solution with fluxes satisfying 



(b) Scalings of the interpolating solution in the 
Einstein Frame. 



64{, 



-)2Ar2 



1. In the left figure, 



c) flow from (^, ^ - 1) ~ (1.22, -0.83) in the IR to (1, -1) in the UV. The right 



'Vor{X5) 

{B,C) - fl™. a™ I 7 1 
figure shows that, in the Einstein frame, the scalings of the (t, x, y)-directions, the w-direction, 
and e't' fiow from f^,^,^) ~ (1.04,0.70,0.70) in the IR to (1,1,0) in the UV. 



33' V33' V33 
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3.3 Interpolating Solutions as Holographic RG Flows 

Here we try to interpret our interpolating solutions via AdS/CFT. To make the argument 
simple we focus on X5 = S^. Consider the standard AdS5/CFT4 for the = 4 super Yang- 
Mills. We can perturb the AA = 4 super Yang- Mills by many relevant operators Oi,(D2,- ■ ■ as 
dS = Jdx^[giOi+g202 + ---]. 

Our solution which flows from the AdS^ x to the scaling D3-D7 solution is dual to an 
anisotropic RG flow triggered by a non-supersymmetric relevant deformation (j2.25p . called Oi. 
This RG-flow eventually ends at the IR fixed point that is dual to the scaling solution. Therefore 
Oi becomes irrelevant in the IR limit. However, our sketch of the holographic RG flow (see Fig. 
[1]) tells us that this flow is unstable and signals the presence of another relevant operator O2 
which becomes relevant even at the IR fixed point. So if we slightly perturb this flow by O2, 
then the RG flow passes near the IR fixed point and eventually goes to infinity. Its asymptotic 
behavior will be derived analytically in the next subsection. However, we can still fine tune such 
that there is no O2 generated (i.e. g2 = 0) to realize the IR fixed point governed by the scaling 
solution. This is exactly what experimentalists usually do to realize an unstable fixed point §| 

A similar situation occurs when we perturb the Heisenberg model by an anisotropic Ising-like 

interaction: H = -JYl{ij){'^^'^^ + '^f'^J + "^f*^!) ~ ^ ^{ij) ^i^h "^^"^^^ D > 0. In this analogy, 
the UV fixed point is the Heisenberg model fixed point and the IR one is the Ising model fixed 
point. This structure of the RG flow is rather generic and is one example of the phenomenon 
called crossover [7J. 



3.4 Exact Solutions without D7-brane Charge 

In our interpolating solution, the flow from the UV AdS^ fixed point reaches the IR D3-D7 
scaling solution fixed point only if it starts along a specific direction. One wonders where the 
flow would end up if its initial direction slightly deviates from the fine-tuned one. A closer look 
at Figure 1(a) tells us that if it deviates from the desired direction slightly to the left, it will 
turn back before reaching the scaling solution fixed point and asymptote to the B < 1 part of 
the black curve given by 9S - 8C - 8 = VSW^T^S (due to the constraint (I3l3]l ). On the other 
hand, if it deviates slightly to the right, it will pass near (but not hitting) the scaling solution 
fixed point and then bend slightly downward and finally asymptote to the B > 1 part of the 
same black curve. 



The curve given by 9B — 8C — 8 = V33i?^ -|- 48 is the solution interpolating between the 
AdS^ in the UV and some other solution in the IR. Since it saturates the inequality (j3.2ip , there 
is no D7 brane charge (/3 = 0): it is a solution of the pure D3-brane system. Not only is it 
important because it gives the asymptotic form of what our interpolating solution decays into 



^It is feasible because as long as the starting point is close enough to the critical line (which corresponds to 
the fine-tuned trajectory that hits the fixed point), the trajectory will stay for a very long time near the fixed 
point to allow the measurements [71. 
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when perturbed by a relevant operator, it is also interesting as a pure D3-brane solution other 
than AdS^. In this subsection, we will solve it analytically. 
First, eliminating C from equations (jS.lip and (j3.12p gives 

^ = 4 + - ^bV^8 + 33S2, (3.25) 

ds A 4 

which can actually be solved analytically. There are two solutions distinguished by the it sign: 



which lead to 



We can choose /iq 
are simply 



c±{s) 



^ and bo 



Co 



C±{s) 



'4 

2 



;8S ^4 



11^ 



aSs 



(3.26) 



+ +^log(e« 



- 1 
4^ + 1 



1 



(3.27) 



Co = 0; then the metric in the Einstein frame and the dilaton 



ds'i 



1 + 1 



+ 



dr^ 



{-df + dx^ +dy^) + 



1 /r^ + l 



dw' 



+ 1 
r4 - 1 



±2 



O0O 



(3.28) 



-N. 



^i*^ = Veto- 

Both solutions in ()3.28p become AdS^ x in the UV fixed point (r — > 00); however when 
going towards the IR, they become singular at r = 1. The solution with " + " (resp. " — ") sign 
covers B > 1 (resp. B < 1) part of the curve - 8C - 8 = V33S2 _^ Even though these 
zero temperature solutions become singular ar r = 1, their corresponding black brane solutions 
at finite temperature are expected to be smooth. An explicit construction of such black brane 
solutions is left as a future problem. Also it would be very interesting to understand what 
anisotropic relevant deformations of the M = 4 Yang-Mills are dual to this background. 



4 Hydrodynamics 

As we obtained the black brane solutions p.27p . we can consider the hydrodynamic behavior of 
their dual field theories from the supergravity side f24]. In this section, we especially determine 
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the shear and bulk viscosities from dispersion relations for the corresponding quasinormal modes 
of fluctuations around the background (j2.27p , following [25l [26] . 

For this purpose, we first reduce part and consider fluctuations around the resulting 
five-dimensional background. The procedure of the reduction and derivation of the linearized 
equations of motion for the fluctuations are summarized in Appendix lA.S.li In these equations, 

(5) 

fluctuation of the resulting five-dimensional metric 5gliJ = H^^, dilaton 5<j) = axion 5x = 
the trace part of the metric vr and the five-form flux 5F^ = /s appear generally. 

In this section, we assume that the w direction is compactified and the neglect the momentum 
in this direction so that we extract the effective 2-1-1 dimensional holographic dual theories from 
the total 3 -|- 1 dimensional ones. In this situation we can choose the momentum in y direction. 
By considering the symmetry of the background, we can decompose the fluctuation into channels 



decoupled from each other in the linear order: 

shear channel Htx,Hxy,Hxr, (4-1) 

scalar channel Hxw, i^-'^) 

sound channel Hu, Hty, Htr, H^x, Hyy, Hyr, H^aiw, Hrr, ^, tt, ftxywr, (4-3) 

"axion channel" Hfw, Hy^, Hrw,rj- (4-4) 



In this paper, we consider the shear and sound channel only, since we can read off the shear and 
bulk viscosity from dispersion relations for the quasinormal modes of these channels. 

Before deriving the dispersion relation, we need to consider the gauge fixing of the fluctua- 
tions. Let us assume that the fluctuations are of the form 

r, = ri{r)e-''^^^-'iy\ vr = ^(r)e-^("*-«^), (4.5) 

Especially, for later convenience, we rewrite some of their components as 

Htt{r) = R^c^Httir), Hab{r) = R^clHab{r), (4.6) 
Hta{r) = R^clHtair), H^^r) = R^ct^^H^Ur), (4.7) 

where a,b = x oi y, and = r'^F{r) and = r^; the function F{r) is defined by (|2.28p . 
As in Appendix IA.8.11 we then take the radial gauge i/^^ = in the following discussion. 
Here we notice that there still exist residual gauge degrees of freedom under an infinitesimal 
diffeomorphism x^' ^ x^' + with = ^M(^)e-iM-w) ; 

H^, ^ H^, - V^e, - V.^M, <f^^- ^^d^'(t>, (4.8) 
V^V- Cpdf'x, vr ^ vr - ^^d^^g^" . (4.9) 
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Here the covariant derivative is defined by the background metric. In [25]! it has been noticed 
that one can derive the dispersion relation for the quasinormal modes by defining gauge invariant 
quantities, instead of performing the gauge fixing completely. Following this, we define the gauge 
invariant quantities as 



Zi(r) = qHUr) + ujH,y{r), (4.10) 



for the shear channel and 



Z^{r) = q'^Huir) + 2qujHty{r) + uJ^Hyy{r) + | (q^^ - ] {H,,{r) + H^^r)), 



(4.11) 



1 

2' 

Z^{r) = TT{r), (4.13) 



Z^ir) = ifir) - -H^Ur), (4.12) 



for the sound channel. Below, we derive equations for these quantities by using linearized equa- 
tions of motion summarized in Appendix[B]and then the dispersion relations for the quasinormal 
modes by imposing an appropriate boundary condition. 

4.1 Shear Viscosity 

From the equations ()B.ip . ()B.2p and ()B.3p . we obtain the equation for Zi as 

14F(r)(q2F(r)-tt)2)-3mVF(r) 121/. (tt,^ - F(r)q^) ^ 
^ 3rF(r)(q2F(r) -tt)2) 36r4(F(r))2 ^ ' ^ - ' 

where the dimensionless frequency tv and momentum q are defined by 
respectively. 

In order to solve this equation, we consider the asymptotic behavior of Zi first. As r — > oo, 
Zi ~ Ar^ + Br~^ with constants A, B generally. For the quasinormal modes, A is set to 
zero, or, in other words, we should impose the Dirichlet boundary condition Zi = at r = oo. 
As r — > //TT, Zi ~ (F(r))^*™/^ and, in order to ensure that only incoming waves exist at the 
horizon, we take the one with negative sign. 

By taking the hydrodynamic limit q ^ 1 and tr ^ 1, we can perturbatively determine Zi as 

Zi = C(F(r))--/2 (l + '-^F{r) + 0(q^ tt))) , (4.16) 

where C is a constant. Here we assumed that q^ and Id are of the same order as usual for the 
shear channel. Then by imposing the Dirichlet condition Zi = at r = cxd, we find It) = — 
Comparing with the hydrodynamic relation for the shear channel 

^ = -^q\ (4.17) 
Ths 
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we can find that the universal bound for the shear viscosity tj to the entropy density s ratio [2^ 
is saturated for the current case: 

T] 1 



s 47r 



4.2 Bulk Viscosity 



From the equations derived in Appendix IB. 21 we obtain the differential equations for Z^, Z^^ 
and as 

,2 ^2^ 



Z:: + ln'(c?c^/3)z; + ^(^-^)z, 



44 11(^3 _ 2^2^,11/3 _^2ri^) 



+ ln'(c?cr^)Z; + ^ f - ^) - = 0, (4.20) 

Z^' + J'{r)Z'^ + g(r)Zo + n{r)Z^ = 0, (4.21) 



where 



625^3 - 3584(g2 _ oj'^ fr'^^ 
" 3r(/i - r"/3)(5;^g2 _ I6rii/3(g2 _ ^2))2 

-20 fi^q'^ {157 - 112uj'^)r'^''/^ + 32^(14^ _ 24u;2)(g2 _ tj2)r-22/3 
+ 3r(^ - riV3)(5^g2 _ I6rii/3(g2 _ ^2))2 ' ^ '^^^ 



_ 3025/1^^"^ - 5^2^2^16/3^333^2 _ i936r-2)(g2 _ ^2^ _ 2304r38/3(g2 _ ^2^3 



2 



9r2(^ _ ^11/3)2(5^^2 _ 16rll/3(g2 _ ^2)) 
3744/xr9g2(g2 _ ^2)2 _^ 5^3^2^5/3(45^4 _ 2541^2^2 _^ ;L93g^2^2) 

+ 9r2(/i - rii/3)2(5^g2 _ I6rii/3(g2 _ ^2))2 ' ^^'^^^ 



,2 , ,2^2 



22(1536ri^(g2 - a;2)3 _ lUOj^r'^'^/^q'^l.q'^ - uj- 



45r2(^ - riV3)(5^g2 _ I6rii/3(g2 _ ^2))2 
22(-^3(425g2 _ 550^^2)^4 _^ 20^^2^2^11/3(83^4 _ 171^2^2 _^ 88^4)^ 

45r2(^ - rii/3)(5^g2 _ I6rii/3(g2 _ ^2))2 ' ^ ' > 

Then, as in the case of the shear channel, all we have to do next is to solve these equations 
by imposing the incoming wave boundary condition at the horizon and the Dirichlet boundary 
condition at infinity to derive the dispersion relation for the quasinormal mode of the sound 
channel. We also take the hydrodynamic limit and assume that q and tr are of the same order, 
as is expected for the sound channel. For = 0, we can see that the nonsingular solution for 
()4.20p is a constant, which turns out to be zero due to the Dirichlet boundary condition = 
at the boundary r = 00. By substituting this into ()4.19p . we obtain = in a similar manner. 
Then from these results and (I4.2ip . by using a similar ansatz Zq = {F{r))~''^^'^YQ(r), we can 
perturbatively determine Zq as 

Zo = C{F{r)r-/' (1 - ^-^^f^^Fir) + 0(q^ m^ qm)) , (4.25) 



19 



where C is a normalization constant. As a result of the Dirichlet boundary condition Zq = at 
r = oo, we obtain the dispersion relation for the sound channel 

Let us recall a hydrodynamic relation for the sound channel in the noncompact (d+l)-dimensional 
spacetime 

V f d-l , ( \ 2 



where Cg and ^ the sound velocity and the bulk viscosity, respectively. Since there are 
two noncompact spatial dimensions for the dual field theory now, d = 2 for the current case. 
Therefore, by comparing (14.26P with (I4.27[) with d = 2, we obtain 

As for (^/rj and Cg of strongly coupled gauge theory plasma in d noncompact spatial dimensions, 
it is conjectured in [28] that they satisfy an inequality 



^>2[-,-4). (4.29) 

In our case, this inequality is saturated. 

5 Entanglement Entropy of D3-D7 Scaling Solutions 

When a quantum system is divided into two subsystems: A and its complement B, the von 
Neumann entropy Sa = — Trp^i log pA (where pA is the reduced density matrix after tracing out 
B) is called the entanglement entropy. The scaling behaviors and certain universal coefficients 
of the entanglement entropy encode important information on the degrees of freedom and non- 
local correlations of the system |29[ [50] . 

For an anisotropic system, an interesting question is "how does the scaling behavior of the 
entanglement entropy depend on the direction along which the subsystems are delineated?" In 
this section, we will study the entanglement entropy of various subsystems A of the (x, y, w) 
space at the boundary (r oo) of the 5D part of the D3-D7 scaling solution (|2.18|) . The field 
theoretical computation of the entanglement entropy is expected to be difficult as the system 
will be strongly coupled. We will instead compute its holographic dual on the gravity side. The 
holographic dual of the entanglement entropy of a subsystem A is given by 

Area'"*'^ 



Here "universal" means the independence from the choice of different dehneations of subsystems. 
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where Area™" is the area of the three-dimensional minimal surface that hves inside the (r, x, y, w) 
space and borders on the boundary dA of the subsystem A [31]. 

After a coordinate transformation r = ^, the metric of the D3-D7 scahng solution ()2.18p 
becomes 

-dt'^ + dx'^ + dy"^ + dz^ dw"^ \ ^ 2 

^2 +^4 



ds^ = [ ^ -^-^ + 7473 ) + ^ ^^i- (5-2) 



For X5 = S^, i?2 = 2\/ttN and ii^ = H^^- We consider the full boundary system given by 



12' 

xe[0,L^], ye[0,Ly], we[0,L^]. (5.3) 

Among the various types of subsystems, we will only study the easiest types: stripes with either 
X or (inequivalently) w restricted to a smaller length. 

5.1 Entanglement Entropy for Subsystem along x-direction 

Let's first consider the subsystem A cut out along the x-direction: 

xG[0,4<iJ, y£[0,Ly], w£[0,L^]. (5.4) 

The three-dimensional minimal surface bordering on dA is given by the embedding function 
z = z(x): 



Area^" = / dx I dy j^^ dw- ^ 



Jo Jo z z 



2/3 



It) + I t^' 1 = R^'LyL^ I dx^Vi + z'\ 







(5.5) 

where d=l-|-l-|-| = |is the total scaling of the boundary system. 

This is a Lagrangian system with £ = R^LyL^j^\/l -\- z''^. The z{x) that minimizes the 
surface area is then given by the equation of motion 



z' = ±\lC-r-l, (5.6) 



z 



2d 



where z^, is the peak of z on the minimal surface, across which z' changes sign. It can be solved 
from £^ = 2 rf* % , which gives z* = % r-^^Y^ ■ The minimal surface is then 

1 /"^ 1 1 

Area™'" = 2i?^L„L^— p- -In with In= / d,,— — (5.7) 

In has a UV divergence at z ^ 0. Imposing the UV cutoff by choosing the lattice spacing 
a for the boundary system, we get 

d-lUa/z^r-^ T{±) i ^ ^ 
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with 



Then plugging the value of and the five-dimensional Newton's constant G5 = 
Gio = Svr^^^, we finally obtain the holographic entanglement entropy for the subsystem divided 
out along the x-direction: 



See- 



V12 



TT 



voKXs; 

71 



" d 



72 



1 



1 



cZ-1 



x5/3 (4)5/3 



(5.9) 



with d = §• 71 and 72 are numerical constants. 

Now let's interpret the result. First, the holographic entanglement entropy is proportional to 
the area of the boundary of the subsystem dA = LyL^ — as expected from the "area law" [29] 
for the entanglement entropy from direct field theory computations. Second, the first term of 
the holographic entanglement entropy diverges and is cutoff-dependent. The scaling of a is given 
by the total scalings of the y and w directions relative to that of the t-direction: (1 + |)/1 = |. 

The second, finite term of the holographic entanglement entropy is more interesting: it 
is cutoff-independent therefore can be compared with the field theoretic computation of the 
entanglement entropy. Its coefficient 72 gives a measure of the total degrees of freedom. The 
scaling of ix is simply the total scalings of the y and w directions relative to that of the x- 
direction: (1 + = |. Since the scaling of the x-direction is the same as that of the 

t-direction, the exponents of a and are the same. 



5.2 Entanglement Entropy for Subsystem along to-direction 

The next easiest subsystem we can consider is to divide along the w-direction: 

we[o,e^<L^]. (5.10) 

The three-dimensional minimal surface bordering on dA is given by the embedding function 
z = z{w): 

r^^ 1 / 

Area's''' = R^LxLy dw^Vz^^ + z'^. (5.11) 
Jo z 

Then we could follow the line of the previous subsystem along the x-direction. The computation 
is straightforward but more complicated so we present here instead a simpler route which utilizes 
the result for the x-direction subsystem. 
The coordinate transformation 



z = z-2, {t,x,y,w) = -{t,x,y,w), 



results in the metric 
ds^ 



-di^ + dx^ + dy^ ^ dz^ ^ dm 



7,2 



+ ^ + ^]+R^ds\^. 



(5.12) 



(5.13) 
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Thus we can simply use the result from the x-direction case, with d = | replaced by dw = 4. 
First, we write down the dictionary between values in the original coordinates and the new one. 



1. In the new coordinates, the full boundary system is 



X G [0, -Lx], 
while the subsystem A is 

X G [0,-Lj, 



y G [0, -Ly] 



y G [0,-Ly], 



w G [0, -Ly,], 



2 2 

w G [0, -(u, < -Ly 



(5.14) 



(5.15) 



2. The lattice spacing in the new coordinates is related to that in the old coordinates by 

2 

a = 03 . (5.16) 

3. The turning point of z is z^, = ^z^. 

Using the result from the x-direction case, we find the minimal surface in the new coordinates 
to be 

1 V^i2h + k] 



Area" 



3-/? LxLy 



1 



1 



(5.17) 



After translated back into the original coordinates, it gives the entanglement entropy of subsys- 
tem along the tf-direction 



See- 



Vol(X5 



N^L^L 



1 



D - 1 



1 



,D-1 



where D — 1 



N'L^Ly 



i2 



(5.18) 



1) 



2 — 2. The negative finite part has the same form as the result for the 
subsystem divided along the x-direction with d = | replaced by =40 

Now let's compare this result with the one for the subsystem along the x-direction. The 
essential features are the same. It is proportional to the boundary area L^Ly. There are one 
cut-off dependent, divergent term and one cut-off independent, finite term. The scaling of the 
cutoff a is given by the total scalings of the x and y directions relative to that of the t-direction: 
(1 + l)/l = 2. The scaling of the £y] is given by the total scalings of the x and y directions 
relative to that of the i(;-direction: (1 -|- l)/(2/3) = 3. Unlike the case for the subsystem along 
the x-direction, since the scaling of the w-direction is different from that of the t-direction, the 
exponents of a and are different. 



^^This might be understood as follows: when all scalings are normalized with respect to the w-direction, then 
the total scaling of the (a;, j/, ■u;)-space is 

A, + A;; + A™ = ^ + I + 1 = 4. (5.19) 
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6 Perturbative Analysis 



In order to know the details of the holographic dual field theories, a basic thing to do is to analyze 
the perturbative spectra around their supergravity solutions. This offers us the information on 
scale dimensions [21 [3]. For example, scalar perturbations are dual to scalar operators Oi in 
the dual field theory. These perturbative modes in supergravity are described by Klein-Gordon 
equations with various masses in the curved spacetime. Since we have the nontrivial dilaton in 
our D3-D7 scaling solutions (j2.18p . it is not clear a priori whether the Klein-Gordon equation 
should be obtained from the string frame metric or the Einstein frame metric. Actually our 
scaling property p.20p of the gravity solutions is only available in the Einstein frame. Also the 
study of the perturbative spectrum is necessary to determine the stability of the background. 
Motivated by these, below we will examine the perturbations around our scaling backgrounds 
(|2J8]) . 

6.1 Description of Perturbations 

Let us analyze the perturbations of bosonic fields around the D3-D7 scaling solutions defined 
by (f238|) . (j2l71) . ([230]) and (plT]) in the Einstein frame . We wih closely follow the analysis of 
AdS^ X in [32l|33j. We will denote the total ten-dimensional coordinates by M, N, .. = 0, 1, 2, •• 
•,9. The five-dimensional Lorentzian spacetime (called M5) is described by the coordinate 
^, z^, ■ • • = 0, 1, 2, 3, 4 and the five-dimensional compact manifold X5 by a, /?,••• = 5, 6, 7, 8, 9. 

As is clear from the IIB supergravity action in our background, the 3-form fluxes F3 and 
H3 are decoupled from the other fields (i.e. the metric, the dilaton, F^, and Fi) thus we can 
concentrate on the latter ones. Then the Lagrangian in the Einstein frame is written as follows 

C = V^(r- \e'^dMXd''x - \dMckd^'ct> - ^FmnpqrF''''^^^^ . (6.1) 

To make analysis more tractable, let us assume = here. Then we can define the scalar, 
vector, traceless symmetric, and antisymmetric spherical harmonics on by , Y^, Y^^^y and 
respectively [32] . Using these spherical harmonics, the metric perturbations Sqmn = hj^N 
can be decomposed as follows [32] 

hi,.)=hl,^Y', h^, = h'Y\ h^^ = B'X, Na/3) = 0'^(i/3)> /^a = vr'^', (6-2) 

where {a(5) denotes the traceless symmetric part. We also denoted all indices of the spherical 
harmonics simply by /. We fix the gauge by requiring 

V"/i(,/3) = V" V = 0- (6-3) 
The perturbations of the dilaton and axion are defined as follows: 

64> = 5x = V^Y^. (6.4) 
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Finally, the perturbation of the 5-form flux = dC^ can be express as follows [32j: 



c — r — r — 

They satisfy the gauge fixing condition V"Ca... = 0, and the self-duality of F5 allows us to 
eliminate 6/,,^2M3M4 and bj^^^^^^. 

Next we substitute (j6.2p , (j6.4p , and (|6.5p into equations of motion of (j6.ip and derive the per- 
turbative differential equations. We omit the details of analysis here and put them in Appendix 
1X1 as many parts of the calculations are essentially the same as those in |32] . 

In the end, we find that the following modes 

Scalar modes : (p^ , (/i^, vr^, 6^), 
Vector modes : (i?^, &^), 

Tensor modes : 6^^, (6-6) 

satisfy free massive field equations which are precisely the same expressions as those in the 
AdS^ X background In the above, the fields in the same parenthesis mix with each other. 
Therefore, these perturbations (j6.6p obey free field equations of motion constructed from the 
Einstein frame metric (j2.18p instead of the string frame metric. 

As an example, let us concentrate on the scalar mode 4>^ . Its equation of motion is written 

as 

+ A) -^'^(1/3)= 0, (6.7) 

where Ox and Dy are the Laplacians of the Lorentzian part M5 and the sphere part S^, respec- 
tively. Using the eigenvalues of ^(q^) , we eventually obtairF^ 

□x - ^^^^^) 0' = 0, (A; = 2,3,4,--.). (6.8) 

On the other hands, the other modes ip^ , r]^ , and h^i, mix with each other and obey equations 
of motion more complicated than those in the AdS^ x case (see Appendix IXj) . 



6.2 Scaling Dimensions and Stability 

We have observed that a large class of supergravity modes (j6.6p . though not all of them, satisfy 
the conventional free field equations with various masses via the Kaluza-Klein compactification 

In other words, for these modes the differences from AdSs x only come from the background metric which 
is employed to write down the free field equations. 

Please distinguish the total angular momentum k of the spherical harmonics of from the number k of 
D7-branes. 
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on . The scalar modes in (|6.6p satisfy the equations of motion of the form (we denote such a 
scalar field by $ here) 

(□^ - m2) $ = 0, (6.9) 

where the Laplacian Ox = g^'^'V^Vi, is constructed from the Einstein frame metric (j2.18p . 
Now let us consider a scalar field <I> on a slightly generalized scaling background 

= ^ (l^ + 72 + ^ j ' (6-10) 

where the scaling exponent u is related to the scaling exponent z in (jl.4p hj u = Especially, 
our scaling background (|2.18p corresponds to u = |. The equation of motion (j6.9p is written as 
follows: 

_ r + ^c^' + + / - ^2 ^ p2 ^2(.-l)^ ^ ^ o_ (g^ll) 

Here uj, p, and pw are the frequency, the momenta in (x, y)- and tf-direction, respectively. 

After redefining the wave function by ^{z) = z^^'if{z), we obtain the Schrodinger-like 
equation 

-^'" + y(z)^' = w^^, (6.12) 

where 

Viz) = ' +p'+pI z^^^-'\ (6.13) 

When z is small, the third term in (j6.13p is small compared to the first term, assuming v > 
Thus, as in the AdS/CFT case (i.e. v = 1), we can expec10 that the stability against the 
(normalizable) perturbations is the same as that of Schrodinger problem with the potential 
= ~^r'- It is well-known that the latter system is stable iff A > 0. 
In this way, we speculate that in the background (j6.10p . the stability condition requires 



^2^2 > _ {^ + ^)\ (g^^4) 
4 

Notice that if we set = 1 in (j6.14p . we obtain rn^W' > —4, which is the well-known Breitenlohner- 
Freedman (BF) bound of AdS^. 

This condition can equally be implied from the behavior of the scalar field near the boundary 
z^O 

(piz) Az^+ + Bz^- + ■ ■ ; (6.15) 

where 



A± = ^± + (6.16) 



We are very grateful to Andreas Karch for illuminating explanations on the stability analysis in Poincare AdS 
spaces. 
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The quantity A± is holographically interpreted as the scahng dimension in the dual anisotropic 
scale invariant theory. Notice that the condition ()6.14p requires that the scaling dimension 
should be real-valued. 

Now let us go back to our D3-D7 scaling solutions (j2.18p . In this case we obtain the stability 
condition by setting u = 2/3 and = j^R^ as follows 

We can apply this condition to the scalar modes in ()6.6[) . As is clear from (16. Sp . all the scalar 
modes cj)^ satisfy this condition. However, we find that one of the infinitely many mixed modes of 
{h^ ,TT^ ,b^) actually has the largest tachyonic mass mpR"^ = —4, which saturates the BF bound 
of AdS^. This occurs only for the second spherical harmonics k = 2. Even though this tachyonic 
mode is stable in the AdS^ spacetime, it seems to become an unstable mode in our D3-D7 scaling 
background as ()6.17p is violated 

Nevertheless, we can replace with an arbitrary Einstein manifolds with the same Ricci 
curvature, keeping the same scaling solution ()2.18p . Define the eigenvalues A of Laplacian of a 
scalar function Y such that —R'^DyY = AY. In this case, if 



A /A 11 , , 

is always satisfied, then the above lowest mass mode becomes stable. In other words, if there is 
no eigenvalue between 

^l^zl^i^ 7.71) < A < ^l±^i^ 16.95), (6.19) 

then the background can be perturbatively stable. Notice that A = 12 saturates the BF bound 
and it is the unstable mode when = S^. It is intriguing whether there exists such a stable 
(Sasaki-)Einstein manifold. 



7 D4-D6 Scaling Solutions 

In type IIA string, the closest analogue to the previous D3-D7 system is the following D4-D6 
system: 



M4xT^ X X^: 


t 


X 


y r 


Wl W2 


Si S2 S3 S4 


N D4: 


X 


X 


X 


X X 




k D6: 


X 


X 


X 




X X X X 



Here {wi,W2) span a two-manifold that supports the D6 flux (we will choose it to be for 
simplicity) and (si, S2, S3, 54) span a four-dimensional Einstein manifold X4 with the same Ricci 

perturbative instability has also been noticed in [2T] for type IIB backgrounds dual to the non-relativistic 
CFT when the 3-form fluxes are vanishing. 
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curvature as that of 5^. However, as we will see, this system does not support a scaling-invariant 
solutions. In this section, we will present a gravity solution that is closest to a scaling solution: 
under (t,x,y,r,wi,W2) — {Xt, Xx, Xy, j, X3wi, X3W2), the line element ds"^ — > X~3ds^. We will 
also give its black brane generalization. 

The fluxes given by these N D4-branes and k D6-branes are: 



(27r)/c 



dwi A dw2 , 



F4 = {27rfN 



1 



Vol(X4) 



-n 



X4 ; 



^2 = 0, 



(7.1) 



where X4 is a unit-radius Einstein 4-manifold (whose Ricci tensor satisfies Rij = Sgij) and 

is its volume-form; L is the periodicity of Wi. This flux profile satisfies the flux equations of 

motion. The corresponding (string frame) metric ansatz is: 



dsi 



(7.2) 



with a possibly non-constant dilaton (f){r). 
For a scaling solution, {A, B,T,Z, 0} are 



A{r) = ailogr + ao, -B(r) = 61 log r + 5o, T = tilogr + to, 
Z = zi\ogr + zo, 4>{r) = ijs log r + 4>o- 



(7.3) 



There are one equation of motion from the dilaton and five from the gravity part. For the scaling 
ansatz, they all reduce to algebraic equations and the solution is easily found to be: 

-3 , , , , , , , 

- log?7s, 



B{r 
T{r 

Z{r 
(t){r 



= (i-5,,)logM4log(3^.*.iv)-l,„g 



= ^%logr + 6o, 



= logr + 1 log (^^«*«) + i log (|.**) . 



6 
2 

3 

^r]s log(r) + ^ log 
77s log r + ^0- 



(7.4) 



3Vol(X4) 



{ris,bo,4>o) are three gauge parameters, rjs corresponds to the gauge freedom of r — r", 60 
corresponds to that of rescaling the {t, x, y) directions, and (j)o gives the string coupling at r = 1 
thus corresponds to rescalings of the r and directions. 
Without loss of generality, we choose 



,, = 2, i^ = il„g(^^^.*Jv)+ilogf, *, = ilog(5^|^) 



3, 34 

+ 2^°^y 

(7.5) 
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The solution in the string frame is 



{-dt^ + dx^ + dy^) + ^ + rt {dwl + dwl) 
r3 



+ Rlr^dsl^, (7.6) 



with Rl = ^Rl = 3YQ||'_y^) • T ■ ^X- Einstein frame, it is 



ds\ = R^ 



1 diT'^ 5 

r3{-dt^ + dx^ + dy^) + — + ri{dwl + d^a) 

r3 



(7.7) 



with R^ = l,R^ = (^vfe)t(i^ • N%fl\ 

This background is no longer scaling invariant. Under the scaling transformation 

T 2 2 

(t, X, y, r, wi,W2) (At, Ax, Ay, -, Xswi, X3W2), 



the metric scales as 



ds^ ^ A 3ds^, 



(7i 



(7.9) 



instead of staying invariant. This is not surprising since the D4-brane theory is not conformal 
in the first place. 

Generalizing to finite temperature, the corresponding black brane solution (in the Einstein 
frame) is 



ds| = R^ 



r 3 {-F{r)df + dx^ + dy^) + 



dr^ 



F[r)r 



with 



F{r) = 1 ir , 
r 3 



- + r3{dwl + dw2) 



+ R^rids^ 



Xai 



where fi is the mass parameter of the black hole. Its Hawking temperature is 

17 3 



The Bekenstein-Hawing entropy is 

• - 



5 = 7 



14 
' 3 



■TJ .iv3.V2-— , 



k 



V3Vol(X4) 

with 7 = 2'3" • 3~3 • 17~3 . TT3 and V2 is the area in the {x,y) directions. 



(7.10) 



(7.11) 



(7.12) 



(7.13) 



8 Conclusions and Discussions 

In this paper, we presented a class of gravity duals of Lifshitz-like fixed points in type IIB su- 
pergravity. They represent backgrounds with intersecting D3 and D7-branes and their Einstein 
frame metrics ()2.18p enjoy anisotropic scale invariance. We also extended them to black brane 



29 



solutions dual to finite temperature theories. Moreover, we showed the existence of solutions 
which interpolate between our anisotropic solutions in the IR and the familiar AdS^ x X5 solu- 
tions in the UV. Then the holography asserts that our Lifshitz-like fixed points can be obtained 
from various four-dimensional CFTs including AT = 4 super Yang-Mills via RG flows. These 
flows are triggered by the relevant and anisotropic perturbation which makes the O-angle (the 
coefficient in front of the topological Yang-Mills coupling F A F) linearly dependent on one of 
the space-like coordinates i.e. 9 oc w. When w is compactified, the perturbation induces the 
Chern-Simons coupling J A A F + ^A^. This theory itself seems an intriguing model worth 
pursuing in a future work, as the equation of motion becomes local in spite of the violation of 
the Lorentz invariance. 

Employing our supergravity solutions we studied the thermal entropy and the entanglement 
entropy to measure the degrees of freedom of the holographic dual theories. We found character- 
istic scaling properties in both quantities. We also holographically computed the shear and bulk 
viscosities. A more general analysis of hydrodynamics with the momentum in the if-direction 
taken into account was left as an interesting future problem. 

Moreover, we performed a perturbative analysis around our solutions and found that a large 
class of scalar modes obey the Klein-Gordon equation in the curved spacetime which has the 
expected scaling property. Also we found an unstable scalar mode when the compact manifold 

is S^. Since this unstable mode occurs only for a 'd-wave' spherical harmonics, might 
decay into a less symmetric Einstein manifold and be stabilized. Thus we have reason to hope 
that there exists a class of (Sasaki-)Einstein manifolds with which our scaling solutions become 
stable. Even the background with = is still useful at least in capturing qualitative 
properties of gravity duals of Lifshitz-like fixed points, with unstable modes simply neglected. 
The construction of manifestly stable and non-dilatonic gravity duals of Lifshitz-like fixed points 
still remains as a very interesting future problem. The three-form fiuxes which we assume to be 
vanishing in our solution might play an important role. 

It is also intriguing to generalize our solutions to other values of p and d in (jl.2p . For 
example, the simplest case p = deserves particular attention. It can be formally obtained from 
our solution (j2.18|) by the double Wick rotation t ^ iw and w — > it. However, the axion field 
X and its fiux become imaginary-valued therefore the solution is not physical in the ordinary 
type IIB supergravity. A slightly better example which realizes the case j? = is a background 
based on D3-D5 systems, where D5-branes are regarded as the baryons [M]. We cannot get any 
consistent solution if we restrict to the ordinary type IIB supergravity because the tadpole for 
the ff-fiux is generated by F-stings which attach to D5-branes and stretch into the boundary 
|34j . To construct a solution in this background we need to add the F-string action as an extra 
term to the type IIB supergravity. Under this slightly unusual assumption, we can indeed find 
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the following black brane solution in the Einstein frame [33] 



2 dp'^ 
' P^A{p) 



g<A(p) = e<^o .6 



(8.1) 



with constant 3-form fluxes H^, and the RR 5-form F^. A derivation of this solution is briefly 
reviewed in Appendix [Cj At zero temperature, this corresponds to the metric (jl.2p with p = 0, 
d = 3 and z = 7. 

Finally, it is also intriguing to apply our backgrounds to realistic condensed matter systems. 
Our D3-D7 model was originally introduced to model the holographic dual of fractional quantum 
Hall effects in string theory [22] (refer to [3^ [371 [38] for other holographic realizations of quantum 
Hall effects). Therefore one of the future problems is to calculate physical quantities such as 
finite temperature corrections to the Hall and longitudinal conductivities in this theory. The 
anisotropic critical points we found in this paper may also be useful to analyze the systems like 
liquid crystals and some anisotropic spin systems. 
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A Perturbative Analysis 



In this appendix, we consider perturbative fluctuations around the D3-D7 background (|2.10p . 
()2.1ip . (|2.17p . and (j2.18p in the Einstein frame with = S^, and analyze the stabihty of 
it by using the hnearized equations of motion. Our analysis below mostly follows \32\ I33j . 
Here we notice that RR 3-form flux and NSNS 3-form flux H3 vanish on this background, 
and their fluctuations do not mix with that of the other fields in the linear order. Thus we 
start with the action (j6.ip . The fluctuation of the metric, dilaton, axion, and RR 5-form flux 
dC^ = along with their decompositions in terms of spherical harmonics on are summarized 
in ()6.2|) , ()6.4p , and ()6.5p . Here we consider the fluctuations satisfying the gauge flxing conditions 
V"/i(q,^) = V^h^a = for the metric and V^CaijKL = for the RR 5-form flux F5 = dC^. For 
simplicity, we denote the background metric, dilaton, axion, and 5-form flux by gMN, (p, X) 
F5, respectively. In this section we deflne a = a/R^ = A/R. 

A.l Some Conventions 

Before writing down linearized equations of motion, we summarize some of our conventions. 

We normalized the e tensor on the flve-dimensional Lorentzian part M5 and the part in 
(I2T8I1 by 

1 



/ 

eoi234 — \/—9M5, e 



01234 



— rP, — ,56789 _ ^ /A i>| 

^56789 — VffS'S' ^ — , ■ (^--LJ 



In this convention, 



We also deflne the ten-dimensional e tensor by 

£0123456789 = 510 = ^01234 • £56789- (A. 3) 

It is also useful to deflne the Laplacian for M5 and by 

a = g'^'V^V,, Dy = 5°^V,V^. (A.4) 

A. 2 Spherical Harmonics on 

In this appendix, we decompose the linearized equations of motion for the fluctuations by using 
the spherical harmonics on S^. Thus we deflne , Y^, Y^^^y and ^[^^p which represent scalar, 
vector, traceless symmetric, and antisymmetric spherical harmonics, respectively. They satisfy 
the transverse conditions 

v-yi = v"y(i^) = v"y[i^] = o. (a.5) 



32 



For the vector spherical harmonics, we define the Hodge-de Rham operator A^, by A^Y^ = 
□j^Y^ — R'^Y^. We can define the Hodge-de Rham operators for the other spherical harmonics 
in the similar manner and the eigenvalues of Ay on the 5^ with the radius R are given b}0 [32j 



A,y^ = n,y^ = -A_Liy/, (fe = 0,1,2,- • •) (a.6) 

A,yi = f - ^ ) = - " ^4^; " (fc = l,2,...) (A.7) 



fc(A: + 4). 

^yYU = (°y - §) = ^ \ Uy (^ = 2,3,4,...) (A.8) 

^vYU = fa. - 41 >^fi«l = (fe = 1, 2, • • •)• (A.9) 



- J Ma/3] " 

A. 3 Five-form Flux Equation 

Let us first consider the perturbation for the RR 5-form flux F^. It satisfies the self-duality 
constraint 

FpQBST = -^^tpQRST^^^^^ FabCDE- (A. 10) 

from which the equation of motion d*F^ = follows automatically. By denoting the fluctuation 
of the 5-form around the background as SF^ = /s, (jA.10|) in the linear order is written as 

f 1^ ABODE f I 1 ^ ABODE TP 

JPQRST - -^(^PQRST J ABODE + ^"^"M ^PQRST J" ABODE 

1 ^ BODE TP UA1A2 /'A ^1^ 

-^^^PQRSTAi FA^BCDEh ^ \ (A- 11) 



By substituting ()6.2p and ()6.5p into this equation and then decomposing it by the spherical 
harmonics on , we obtain five equations: 



= 0, (A.12) 



[bU,.,.,. + ^,.,.,.,yy] = 0, (A.14) 

1 

'T? Vl , v\i'2ul , 717273V7 T^/ _ n (' A 1 ^^^1 

^/*l>2/<3^[a/3] Y2 '^i^2/'3 %i^2^a/3 ^ 7i [7273] ~ (^A.iOJ 

[bl,2,-s + ^,.,2,7'"'^M = 0. (A.16) 



Here a = a/R^ = 4/R. Now, we can simply solve (IA.14P and ()A.16p algebraically, assuming 
k > 1. Then, we obtain the following three equations from ()A.12p and ()A.14p . (jA.lSp and ()A.16p . 



Please distinguish the total angular momentum k of the spherical harmonics of from the number k of 
D7-branes. 



33 



and (|A.15p respectively: 



= 0, 



2 2 

p^bj, - v^v.bi + A,5^ - di?^] yi = 0, 



3V„„ bi 



f 1^2 hi 



a/3] 



0. 



(A.17) 
(A.18) 
(A.19) 



A. 4 Einstein Equations 



Let us next consider the perturbation for the Einstein equation. The Einstein equation itself 
can be obtained from the action ()6.ip as 



Run - ^9MN {R - ^e^^Opx^^x - ^9p</'5^</' - -j^FpQnsTF'''^''''A 



1 

T 

\e^'^dMXdNX - 7;dM4>dN(t> 



4-5! 



4-4! 



'MPQRS-f N 



PQRS 



0. 



By using 



and 



(A.20) 

(A.21) 
(A.22) 



derived from the trace part of ()A.20p and the self-duality condition (jA.lOp respectively, we can 
reduce the Einstein equation ()A.20p to a simpler form: 



1 



1 1 

Rmn - -e'^'^'dMXdNX - ■^dM(t>dN(t> ^ ^, 



mpqrsFn 



PQRS 



0. 



(A.23) 



By linearly perturbing ()A.23P and decomposing it in terms of the spherical harmonics on S^, we 
obtain some equations for the fluctuations around the background. We summarize the resulting 
equations below. 



A. 4.1 a/5 Components 

From (a, /5)-components of ()A.23p . we obtain the following four equations: 



□x + - — 



2 

R\ 

v^i?^V(„y/^ = 0, 



V(„v^)y^ = 0, 



5 

1 ^ f 4 r 1 , r„ d, ^ d^ r 

10 25 ^10 2 ^ 5 



(A.24) 
(A.25) 
(A.26) 

(A.27) 



34 



A. 4. 2 fia Components 

In a similar way, from a)-component of ()A.23p . we obtain the following equations: 

Y^ = 0. (A.28) 



P2P3P4 



1 



1 



a 



UfJ, 



96 

v„y^ = 0. 



P1P2P3P4 



(A.29) 



A. 4. 3 /iz^ Components 

Finally, from (/i, i/)-component of ()A.23p . we obtain the following equations: 



1 



~ 2 

-i^d.c^d.^ - ^-d^d,^ - ^5,,e^iP2^3P4P5v^^6^,,3,,,3y^ - ^g,.h + V- 



T 



(A.30) 



A. 5 Dilaton and Axion Equations 

Let us then move to the linear perturbation for the dilaton and axion equations of motion. They 
are given by 

5m(^/^/''^9^0) - ^/^e^^/^^^MX^ivX = 0, (A.31) 
9M(^/^e2'^/^^9^x) = 0. (A.32) 

respectively. When fluctuating around our background, these equations are rewritten as 
+ Uy)^ + ]^g'-dr(pdr{h + vr) - {dr(^){Vy) - {Vidjct>)h'J 

-2(^e2 + e^'t'h^^id^xf - 2e^'t> g^^ {d^xW^r,) = 0, (A.33) 
+ ay)rj + 2idr(P){drr])g'' + 2g^^ idy,^)id^x) + + vr)a^x 

-(V^/i'^"')(a^x) - h'^iVidjx) - 2W-^{dr4)[d.^x) = O. (A.34) 



A. 6 Spectrum for Decoupled Modes 



As we derived the linearized equations of motion for the fluctuations in the appendix IA.31 IA.41 
and IA.51 we then analyze the spectrum for them. We start with the one for those modes which 
do not mix with other modes in a complicated way. The analysis turns out to be essentially the 
same as the case of AdS^ x [32j . 
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A. 6.1 4)^ Mode from 



First we consider the scalar mode /i(Q^) = (j) ^(q,/3) from the fluctuation of the metric. It obeys 
the equation of motion (|A.24p and, by using (|A.8p . we obtain 

(^□, - / = 0. (A; = 2,3,4,...) (A.35) 

Thus we find its mass m? = which obviously satisfies the stability condition. 

A. 6. 2 Mode from Cf^ua/B 

For bjj^j^ from the fluctuation of C^,yai3 > (|A.19p leads to 

(Max, + A,)6;,,y[^^]=0. (A.36) 



The Maxwell operator Max^; on M5 is defined by Max^.^^ = Dx^^^ — V'^V^bj^ for the vector 
6^ and we can generalize the definition for teonsors. Since the mass for this mode is given by 



m? = —Ay = ^^^2 = 1, 2, • • •), this mode turns out to be stable. 



A. 6. 3 6^ and B^^ Modes from g^a and C^aP'^ 

Let us next consider 6^ and from the fluctuation of g^a and C^aP'y respectively. From (jA.lSp . 
()A.28p . and ()A.16p . we obtain the equations for these modes as 

((Max, + Aj,)6^ - aBl)Yi = 0, (A.37) 
((Max, + I\y)Bl - ^{Max - \)h'^)Yl = 0. (A.38) 

or, by denoting in a different expression, as 

The eigenvalues of 2 x 2 matrix in the second term are 



m| = A,-^±V^-«2A, . (A.40) 



More explicitly, by using ()A.7p . we get the masses for these vector fields 

2 9 (A; + 3)(/c + 5) , 
ml = ^^, ml = ^ A (fc = 1,2,3,- ••). (A.41) 

Therefore we can find that these modes are stable, too. 
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A. 7 Spectrum for Mixed Modes: , tt^, and 

Let us determine the spectrum for the scalar perturbation , vr^, and here. We first assume 
/c > 2 and then find = -^n^ from (1X26]) . Thus we can rewrite ()A.17p and (|A.27p as fohows: 

The matrix in the second term is diagonalized and the eigenvalues are given by 



mi = Uy- ± Ja'^ - ^a^Uy, (A.43) 



or, more explicitly, by 



ml = — (A.44) 

ml = i^±^l^. . (A.45) 

Even for /c = 0, 1, we can see that the expression ()A.45P is correct. Thus we find the lowest mass 
in this mode is = —-^ when k = 2. This violates the stability bound. 

A. 8 Spectrum for Other Modes: /i^jj,, 99, and rj 

Now, the remaining modes are /i^^,, (/?, and 77. Since the analysis of massive modes looks highly 
complicated, below we would like to consider only zero modes on of h^u and ^p, tj. They 
are useful to the calculations of the viscosity in section HI For this reason, we generalize the 
background to the black brane metric ()2.27p and write down the linearized equations for the 
fluctuation around it. By taking — > 0, we can reproduce those for the background metric 
(f2T8]) . 

A.8.1 Zero Modes on 5^ 

We concentrate on the zero modes on i.e. A; = modes of the spherical harmonics. Notice 
that in this case, we have 6 = and 6/ji/x2/i3At4 is expressed in terms of h and vr as 

2 

5V^j^6^2Ai3M4M5 ~ ^^A'1M2A»3M4A»5 ~ 1 (A. 46) 

from (jA.lOp . From (|A.27p . the mode vr satisfies 

32 \ 

□x - ;^ j ^ = 0. (A.47) 
It is also useful to define the Weyl shifted metric Hfj,i, by 

Hfj,u = h^j,y + ^g^.v'K, (A.48) 
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which corresponds to the metric perturbation around the five-dimensional background obtained 
by the reduction of S^. Then the Einstein equation ()A.30p can be rewritten as follows: 



4 e^'f' 1 

+ j^H^, - e^Hd^xduX)'P - ^(^^.X^.r? + d.xd^r,) - -{d^<t>duV + d^cj^d^^) = 0. (A.49) 

Here 5R^I^} is the perturbation of the purely five dimensional Ricci tensor (neglecting the 
contributions) due to the metric perturbation H^jj. Notice that when ^ w, r, we obtain the 
simple Einstein equation 6Rliu + -^H^^ = 0. This simplification is applicable, for example, to 
the shear modes Htx and Hxy 

For the dilaton equation of motion (jA.SSP with /c = 0, by using the results 

2 2/-, ^^ \ f, , 5/i 



22/1 ^^ 

= (1 - ^ 

and the Weyl shifted metric H^^, we can rewrite it as follows: 

□xV^ + ^g^'idrH) - ftH'' - + Ryy) - UH^^ - frH^' 

+(5'^e^'f'{H'"'" - 2g'"'"^ - -irg'"'") - 2(3e^^ g'"'" d^i^ = 0, (A.51) 

3 

where H = H'^ = h + ^TT. 

On the other hand, for the axion equation of motion ()A.34j) with /c = 0, by using the values 

2k 

Vrd^X = ^wdrX = -17^, (A.52) 
SLr 

it can be rewritten as follows 

We can take the radial gauge H^^ = and have 5 physical modes for H^^. Then the dilaton 
and axion equations of motion ()A.5ip and ()A.53P become 

□x<^ + ^g'^'^drH) - ftH'' - + Ryy) 

6r 

f ttww I ^ ^/ tt f , r, XX f I WW j: I rr f \ 

-jwH +3^(5 + Jx + g jw + g jr) 

^p2^2<P^jjww _ 2ga>w^ _ ^T^g^n,) _ 2^6^''' g"^"" 8^7] = 0, (A.54) 
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and 

OxT] + Y/"drV + fig'^'^d^ (^ip + f ) - /3V^^^" = 0. (A.55) 

In summary, we need to solve the zero mode equations of motion ()A.46p . (|A.47p . ()A.49p . 
()A.54p . and ()A.55P to find the variables h^yp^j, it, Lp, rj, and H^y. 



B Linearized Equations for the Shear and Sound Channel 

In this appendix we summarize the linearized equations for the shear and sound channel. They 
are useful to derive the differential equations for the gauge invariant combinations in section [H 
For this purpose, we consider the equations (|A.46p . (pV.47p . ([X!49i) . ()A.54p . and (|A.55p derived 
in the appendix lA. 8. II Then we substitute the fluctuation corresponding to the shear channel 
(j4.ip or sound channel (j4.3p into them and then derive explicit expressions for the linearized 
equations of motion. We summarize the resulting equations below. 

B.l Shear Channel 

By considering the fluctuations Htx and H^y which correspond to the shear channel in the radial 
gauge Hfj^r = 0, we obtain the equations for them 

H'l, + \x^{cl^l^)H[, - -^{qHtx + ujH^y) = 0, (B.l) 

H'^y + \^{4cl'^)Ky + -^{qHtx + u^Hxy) = 0, (B.2) 

qH',y+u:^H[, = {), (B.3) 

from (x,y)-, and (j;, r)-component of (|A.49p . respectively. We can directly show that two 

of the three equations are independent. 



B.2 Sound Channel 

Let us next consider the sound mode. The corresponding fluctuations are Htt,Hty, H^x, Hyy, Hww, 
and //xi/x2M3M4A«5 ™ radial gauge. In this case the linearized equations (|A.47p . (|A.49p . and 
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()A.54p are summarized as 



H'{, + In' (c?4/3) - In' (q) {H' + H[, 



yy 

H'l, + In' (cfc^/s) //^^ - ln'(c,) (//^ - //^^ - H^) + ^{^-^^ H^, = 0, (B.5) 

H';^y + \^{44!^)H'yy - ln'(c,.)(^; - H'y.y ■ H^^ 

H — 2 \~^yy + ~^^ty + -^{Htt — -f^ii) ) = 0; (B-6) 
V Cj c^. / 

<^ + ln'(c?c8/3)//;^ - ln'(c2/3)(i/; - H'yy - H^^ 

+ ln'(c?c-2)//;, - In' (^c^' - i/^;^) = 0, (B.8) 
H'{y + ln'(ci^/3)i/;^ + = 0, (B.9) 



{H'yy + F;,) + ^H[, + ln'(c,c,-^)(i/,, + 



+^ ln'{clc^^)Hty + ln(c2/3) (^^ - ^H^J^ = 0, (B.IO) 
{Hi, - H[,) + -4^L - ln'(c,.qi)Fii - ln'(c2/3) L - i//^^") = 0, (B.ll) 
^" + ln'{cyj')^' - ln'{cl/'){Hl, - H'yy - H[^ + ^ f ^ - 4) 

"9^ Y ~ 2"""") + 27r28/3(rii/3-^) " °' 
vr" + ln'(c?4/3)vr' + ^ (i^ - vr - ^tt = 0, (B.12) 



where i^ii = i?^^. + H^^. 



C D3-D5 Scaling Solution with F-string Sources 

Here we briefly review the D3-D5 scahng solution in type IIB supergravity with F-string sources 
[35j . This solution can be regarded as a back-reacted solution dual to a homogeneous baryon 
condensation in = 4 super Yang-Mills. Each baryon (= D5-brane) situated near the horizon 
r = carries N F-strings which extend into the boundary r = 00 due to the string creation |34j . 
Below, we work in the Einstein frame and the supergravity action is obtained by rewriting the 
one in the string frame ()2.ip by the Einstein frame metric Gf^j^ = e~'^^'^G^j^}^^ . 
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For this brane setup, it is appropriate to use the fohowing ansatz for the metric 

ds'^ = -e^"-^ A{r)dt^ + e*^'^) {^{dxif^ + e^^'XA^^ {r)dr'^ + r'^ds\^), (C.l) 

where represents a Einstein manifold with the same Ricci curvature as the unit radius . 
For fluxes, the ansatz is 

F3 = et'''^'"^ • h{r) ■ dxi A dx2 A dxs, 

F5 = r5e^M/(r)(17x,+*J^Xs), (C.2) 

where is the volume form on X5. We also assume that the dilaton is dependent on the 
radial coordinate r only. If we treat this solution within supergravity, we will have a tadpole 
of the NSNS 3-form flux H-^ as is clear from the equation of motion in the presence of the 
Chern-Simons term 

d{e-'t'*H^) = F^^F^, (C.3) 

where and F5 are sourced by the D5 and D3-branes, respectively. Thus we cannot construct 
consistent supergravity solutions under the assumption of symmetry of spacetime implied by 

(|aTl) . 

To resolve this problem, we notice that the effective F-string charges are generated from the 
D3-brane and D5-brane charges via the Chern-Simons term 

Scs = -^ J CiAHsAFs. (C.4) 

Therefore, to resolve this problem, we treat the created F-strings as the probe action and add it 
to the supergravity action. The probe action for a single F-string in the Einstein frame is given 
by 

Sstr^ng = / ^rdae^ y^G^ + ^ / dTdaBf,^^rX^'^^X^'. (C.5) 

By identifying t = t and a = r and taking the sum over infinitely many F-strings, we obtain 
the probe action for the created F-strings 

Sstring(i) = j dxidx2dX'i^x^^ P Sgtring, (C.6) 

i 

where i labels i-th F-string and the density of F-strings p is assumed to be constant. In the 
following, we solve the equations of motion derived from the supergravity action with this probe 
action added. 
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C.l Equation of Motions 

The Bianchi identity for and the equation of motion for are written as 

9^(r5ei'=W/(r)) = 0, a^(ei^W/i(r)) = 0. (C.7) 

Here we notice that the equation of motion for F3 is automatically satisfied for the ansatz 
introduced above. From these equations, we can define the constants F and H as follows: 

/(r) = F-r-^e-i'^^'^), /i(r) = • e-f^'W. (C.8) 

Then, from ()C.3p . the density of F-strings turns out to be 



p=— ^. (C.9) 

To derive the dilaton equation of motion, we notice that the F-string action now looks like 

Sp-string = j dtdre^^-GE = j dtdre^+^+^ . (C.IO) 

Thus the equation of motion becomes 

a.(0'(r)A(r)e2'=M+i''M+^r5) = ^e^+^+^ + ir5e<^M+i^M+3^M+^/i(r)2. (C.ll) 

We can also derive the Einstein equations for the type IIB supergravity with the F-string 
action. Combined with (jC.lip . we can summarize the equations of motion as follows: 

[^(log^ + uyr^e-2'+^'+-2^]' = ^FHe'2'^+"2-+h + 1^^-56-2^+1^+^ + ^H^^e^+^'-l'^'^- , 

(C.12) 

L J 8 2 4 ' V • ; 

[(c + 21ogr)'Ar5et''+2^+^]' = -1^^65-^+5^+^ - iFV-^e-^'^+f ^+^ 

8 2 

+li/V5e^+3^-i''+5« + 8r3e2-+i^+l", (C.14) 

[(log A + u-b- 20)'^r5ei^+2^+3"]' = 0, (C.15) 

2{(P'f + 3(6')^ + 10- - 3b'c' + 66" + 10c" - 3b'u' - 5c'u' - 10- = 0. (C.16) 



It is also useful to derive the following equation from a linear combination of ()C.12p , (jC.lSP , 
and (IClil) 

[(log^ + n + 6 + 2c + 41ogr)'^r5ei*+2^+^"]' = IGr^e^^+t^+l". (C.17) 

In the above discussion, we have derived five equation of motion for five variables A{r), 
u{r), b{r), c(r), <j){r). However, we can eliminate one of them, say u(r), by the diffeomorphism 
r ^ r = f{r). In this sense, the independent degree of freedom under the symmetry ansatz is 
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the four variables. Thus we should show the five equations of motion are not over constrained. 
Indeed we can show the following identity from the four equations of motion ()C.12p . (jC.lSp . 
(IClill . and ([Cl5|) 

,10^(,)2^36(.)+4cM+nM . E,,^,,ra^ntir)]' = 0, (C.18) 

where Econstrainti''") is the left-hand side of (|C.16p . This guarantees that the fifth equation of 
motion Econstraint{r) = is satisfied if it is vanishing at any particular value of r. 

C.2 Scaling Solutions 

Let us assume the following simple scaling profile for the unknown functions: 

u{r) = ui log r + uq, 6(r) = 6i log r + 6o; c(r) = ci log r + cq, i;^(r) = logr + ^q. (C.19) 

At first we also assume A{r) = 1, which corresponds to the extremal case. 

Using the three equations of motion (jC.14p . (jC.15p . and (|C.16p . we find that the following 
coefficients ui,bi,ci and (pi satisfy equations of motion 

^1 = 7^1, ^i = yf, ci = -2, '/'i=3y|. (C.20) 

By substituting (IClOll with (10201) into (10121) . (IOT3|) and (lOTil) . we can find four solutions 
for the pair {F, H) in terms of uq, bQ, cq and (pQ as follows: 



F = ±e 



H = ^^~^ei''o-i-°-^^°-F, (021) 
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and 

F = ±e 



\/83 + ^/l081 
V3 ■ 



^{I081+5.^|.o-|co-i*o.^. (C.22) 

A A ^ ■' 



" 44 

Thus, in order to realize the fluxes with arbitrary values and sign, we have only to choose one 
of the four solutions and tune cq and (po appropriately. We set mq = = below. 

As we constructed the extremal solution, generalization to a black brane solution is straight- 
forward by considering more general functions A[r). Let us again assume the ansatz ()C.19p with 
the values ()O20p . Then it immediately follows that all equation of motion are satisfied only if 

A{r) = 1 - (C.23) 

where M is an arbitrary constant related to the ADM mass of the black brane. We also notice 
that the other profiles are the same as the extremal solutions. 
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In summary we obtain the following metric in the Einstein frame and the dilaton 

ds'^ = -r''M(r)dt2 + ^feirff 2 ^ + {dn^f^ , e'^(^)-'^« = r'^' . (C.24) 



and, after the redefinition of radial coordinate p = , we reach the expression ()8.ip . 
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